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ABSTRACT 

Numerical hydrodynamical simulations have proven a successful means of reproducing 
many of the statistical properties of the Lyman-Alpha forest as measured in high 
redshift quasar spectra. Pseudo-hydrodynamical methods based only on simulating 
the dark matter component have been claimed to yield a comparable level of success. 
We investigate the degree to which two pseudo-methods, with and without allowing 
for a pseudo-gas pressure, are able to match the predictions of fully hydrodynamical 
plus dark matter simulations. We also address the requirements for convergence to the 
statistics of the spectra and the inferred properties of the Lyman-Alpha forest as a 
function of resolution and box size. Generally we find it is possible to reach agreement 
with full hydrodynamic simulations at the 10% level in the cumulative distributions 
of the flux and absorption line parameter statistics for readily achievable particle and 
grid numbers, but difhcult to do much better. 
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1 INTRODUCTION 



Numerical simulations of structure formation in the universe 
incorporating hydrodynamics in Cold Dark Matter (CDM) 
dominated cosmologies have proven very successful in repro- 
ducing the statistical properties of the Lya forest as mea- 
sured in high redshift Quasi-Stellar Object (QSO) spectra 
(Cen et al. 1994; Zhang, Anninos & Norman 1995; Hern- 
quist et al. 1996; Zhang et al. 1997; Bond & Wadsley 
1997; Theuns, Leonard & Efstathiou 1998a). The high level 
of statistical agreement suggests that the models are cap- 
turing the essential physical nature of the absorbing struc- 
tures (Meiksin et al. 2000). Much of the structure of the 
Intergalactic Medium (IGM) may be understood as a con- 
sequence of the spatial coherence in the statistical proper- 
ties of the density fluctuations of the dark matter alone in 
a "cosmic web", as emphasized by Bond & Wadsley and 
supported by the slow comoving evolution of the integalac- 
tic filaments found by Zhang et al. (1998). Gnedin & Hui 
(1998) claim that in fact it is unnecessary to incorporate full 
hydrodynamics in the simulations to reproduce the proper- 
ties of the LyQ forest. Gnedin & Hui introduce instead a 
pseudo-hydrodynamical scheme into an N-body calculation, 
reporting comparable results in the distributions of flux and 
H I column density to full hydrodynamical -|- N-body simu- 
lations for a standard CDM model. Pseudo-hydrodynamical 
methods based on PM alone have been used by Petitjean, 
Miicket & Kates (1995) and Croft et al. (1998). 

Because of the large computational overhead in incor- 
porating proper hydrodynamics into an N-body code, con- 



siderable beneflt would be obtained if the hydrodynamics 
computations could be avoided. It would then be possible to 
perform a wide range of cosmological models at only mod- 
erate computational expense, as well as test systematics on 
the properties of the forest: the effect of box size, resolution, 
and even physical effects like variations in the character of 
the background UV radiation field that ionizes the gas. It 
would be of additional theoretical interest to infer the degree 
to which the properties of the forest are hydrodynamical in 
origin, as opposed to being solely a consequence of the clus- 
tering of dark matter. 

To these ends, we perform simulations of the Lya forest 
for a variety of cosmological models using two approaches: 
a pure gravitational computation using a Particle-Mesh N- 
body code (PM simulations) , and a pseudo-hydrodynamical 
technique, based on the method of Gnedin & Hui, incor- 
porating an effective gas pressure force computed from the 
local dark matter density (HPM simulations). 



2 THE MODELS AND SIMULATIONS 

We consider two sets of models. The first is based on sim- 
ulations in a flat CDM model with a cosmological constant 
(ACDM ) using both the PM and HPM codes. These are per- 
formed in order to test convergence of the results for differing 
box sizes, particle numbers, and grid resolutions. A second 
set is based on 4 CDM models in Machacek et al. (2000), 
who performed computations using full hydrodynamics. We 
duplicate their models in order to compare our resulting 
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Figure 1. Non-dimensional linear power spectra A'^(fc) = 
P(k) / (2iT^), where P{k) is the power spectrum, for the four 
cosmological models examined. Shown at z = 3. 



positions by assigning particles to a regular cartesian grid 
using the C loud-In-Cell (CIC) scheme (Hockney & East- 
wood ( 1988| )). The Fourier Transform is taken and the force 
computed using the kernel k/k^. We use a second order leap 
frog method to integrate the equations. The relevant posi- 
tions are predicted at a half time step and used to calculate 
the accelerations which modify the velocities. The time step 
is dynamically chosen as a small fraction of the free-fall time 
with a maximum size of A log a = 3%. 

The initial conditions were created by displacing the 
particles from a uniform grid using the Zel'dovich approxi- 
mation. 

Our implementation of the HPM code differs slightly 
from that in Gnedin & Hui. Firstly we use 2 sets of par- 
ticles, one of which feels only the gravitational force (the 
'dark matter') and one of which additionally feels a pseudo- 
pressure force (the 'gas'). The gravity is computed as before. 
To compute the 'hydro' forces we bin the 'gas' particles onto 
a grid using CIC as for the gravity. We compute the tem- 
perature using the equation of state 



r = To I ^ 



7-1 



(2) 



where To and 7 are obtained from Hui & Gnedin (1997). 
The enthalpy per particle 



7 



(3) 



statistics for the absorption features with theirs. The mod- 
els investigated, in additional to the ACDM model, are: a 
standard critical-density flat CDM model without a cosmo- 
logical constant (SCDM), an open CDM model (OCDM), 
and the standard CDM model but with the power spectrum 
of the density perturbations tilted (tCDM) to match the 
normalization on large scales as determined from the COBE 
measurements of the Cosmic Microwave Background (Bunn 
& White 1997). The model parameters are shown in Table 
|l| As in Machacek et al. , the initial data were generated 
using the BBKS transfer function (Bardeen et al. 1986) 
to compute the starting redshifts and the unconstrained ini- 
tial particle positions and velocity perturbations appropriate 
for each model. The non-dimensional linear power spectra 
A^(fe) are shown at 2 = 3 in Figure |l|. 

To account for the degree of ionization of the hydrogen, 
it is also necessary to include a background UV radiation 
field. We essentially circumvent this by fixing the mean op- 
tical depth, as defined by -fa = — log < flux >, at fa — 0.30 
at z = 3, as in Machacek et al. 

The evolution is computed using the PM code de- 
scribed in de tail in Meiksin, White & Peacock (1999) and 
White ( [l999| ). The computational volume is chosen to be 
periodic with side length 1, and we take as our time coordi- 
nate the log of the scale factor, logo, where a = (1 + z)~^. 
Velocities are measured in units of the expansion velocity 
across the box aHLhox, where H is the Hubble parameter 
fixed by the Friedmann equation 



= H, 



0,1 



(1) 



in terms of the Hubble constant Ho and the densities in 
matter, cosmological constant and curvature, in units of 
the critical density. The density is defined from the particle 



is then computed, Fourier transformed and smoothed by a 
gaussian of half a grid cell (to supress high frequency noise) . 
The 'hydro' force is then calculated from the inverse Fourier 
transform of ikip and added to the gravitational force for the 
'gas' particles. 

In practice, To will depend on the reionization history of 
the IGM. In particular, if He 11 is reionized (to He 111 ) late, 
lower values of To will result prior to He 11 reionization and 
higher values in rarified gas after (Meiksin 1994). The effect 
of reionization on the thermal properties of the gas will only 
be known once the nature of the ionizing sources has been 
established, and may eventually require the incorporation of 
radiative transfer into the simulations. 

We consider a variety of particle number and grid zone 
combinations, as well as several box sizes. The PM simu- 
lations are summarised in Table ^ For the HPM simula- 
tions, two combinations are used: {Np, Ng) = (128, 256) and 
(256, 512). For HPM the dark matter and the gas are rep- 
resented by an identical number of particles, each equal to 

Given a set of final particle positions and velocities, 
we compute the spectra as follows. First the density and 
density-weighted line-of-sight velocity are computed on a 
grid (using CIC interpolation as above) and Gaussian 
smoothed using FFT techniques. This forms the fundamen- 
tal data set. A grid of sightlines is drawn through the box, 
parallel to the box sides and along each sightline a new ID 
grid of density and velocity are obtained from the funda- 
mental grid using CIC interpolation. The position in veloc- 
ity units and the temperature are computed on this ID grid 
using the cosmological model and equation (2) for the equa- 
tion of state. 

Using this ID grid we then integrate in real space to 
find t{u) at a given velocity u. Specifically we define 
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where mo = xaHLho^ + uios and b = •y/2A;sT/mH is the 
Doppler parameter, where mn is the mass of a hydrogen 
atom. The flux at velocity u is exp[— r]. The integration vari- 
able x indicates the distance along the box in terms of the ex- 
pansion velocity across the box. The constant A is then itera- 
tively adjusted to obtain a predefined fa = — log {exp(— t)). 

We generate sample spectra from the simulated data 
along random lines-of-sight through the computational 
boxes over the length of each box. In order to make a fair 
comparison with the data in MOO, both sets of spectra are 
rebinned to a constant pixel resolution of A/ A A = 74000 
and gaussian smoothed to mimic a spectral resolution of 
8kms~^. These parameters were chosen to approximate 
what may be achieved using the Keck HIRES. 

There are several measurements of the Lya forest which 
may be used as a basis for comparison. The most fundamen- 
tal is simply the distribution of flux per pixel in the syn- 
thesized spectra. The fluctuations in the flux on different 
velocity scales provide additional constraints, which may be 
related both to the density fluctuation power spectrum and 
to the thermal widths of the lines. We quantify the fluctua- 
tions in three ways. We compute the power spectrum Ppik) 
of the flux, an analog of the density power spectrum. (For 
the Ppik) estimates we use the raw spectra from the PM 
and HPM runs.) The widths of the spectral features may 
be directly characterised using wavelets, which quantify the 
changes in a given spectrum on being smoothed from one 
velocity resolution to another (Meiksin 2000). A multiscale 
analysis based on the Daubechies wavelets effectively per- 
forms a weighted smoothing of the spectrum over successive 
doublings of the pixel width. We use a Daubechies wavelet 
of order 20. Lastly, the most traditional method to quan- 
tify the flux fluctuations is by decomposing each spectrum 
into a set of Voigt line profiles. We utilise AutoVP (Dave et 
al. 1997) to perform such an analysis. Each sample spectrum 
produces an average of on the order of 900 lines per unit red- 
shift. Typically (3 — 6) x 10'* lines are analysed per model. 
The absorption lines are characterised by their H i column 
densities and Doppler parameters. The H i column densities 
reflect both the density and physical thickness of the struc- 
tures which give rise to the absorption features. The Doppler 
parameters reflect the degree of line broadening due both to 
the thermal motion of atoms and any internal peculiar ve- 
locities within the structures, and may provide probes of 
the temperature and velocity history of the universe at high 
redshifts. 

To ensure our results are not affected by sample vari- 
ance, we performed two independent simulations of some of 
the models. We found convergence on the statistical prop- 
erties could be achieved with a single realization provided 
the statistics were averaged over a sufficient number of lines- 
of-sight. We found 3072 line-of-sight (each one box side in 
length) was adequate for the PM simulations, and 1024 for 
HPM. 



3 CONVERGENCE TESTS 

Simulations of the Lya forest are restricted by two size con- 
straints. High spatial resolution is required to resolve the 
structure of the absorbers, while a large box size is required 
both to capture the large-scale power which will affect the 
velocity widths of the absorption lines as well as to obtain 
a fair sample of the universe. Although a small box may be 
adequate to resolve the line structure, if the box is too small 
the fluctuations corresponding to the scale of the box will 
become nonlinear, in which case the results are no longer 
representative of the cosmological model simulated. 

Resolution effects were previously studied using full 
hydrodynamics by Theuns et al. (1998b) and Bryan 
et al. (1999). Using P^M-SPH, Theuns et al. found for 
SCDM that a (comoving) box size of at least 5.5 Mpc using 
64'' SPH particles and 64^ dark matter particles is required 
to converge on fa and the Doppler parameter distributions 
at z = 3. This corresponds to a mean proper interparticle 
separation of 21 kpc, or about half this in structures with an 
overdensity of 10, typical of the filaments with H i column 
densities of lO^"* - 10^^ cm~^ (Zhang et al. 1998). Bryan 
et al. found that a box size of at least 9.6 Mpc was needed 
to reach convergence on the Doppler parameter distribution 
at 2 = 3 using a PM-PPM code (Kronus) with 128^ grid 
cells, corresponding to a proper resolution of 19 kpc. 

Neither group was able to maintain a resolution of 
10 kpc while increasing box sizes beyond 5 Mpc due to 
memory limitations, so that the contribution of larger scale 
power at high resolution was left unexplored. We also note 
that SCDM has less large-scale power than more viable cos- 
mological models like ACDM , so that the resolution limita- 
tions are least restrictive for SCDM. It is possible to achieve 
much higher resolution using PM or HPM. In this section 
we present the results of a series of tests of the convergence 
of the PM and HPM simulations for ACDM as the box size, 
particle number, and number of grid zones are varied. The 
simulation results are all shown at 2 = 3. 



3.1 Flux distribution 

A comparison of the cumulative distributions of flux per 
pixel for several resolutions and box sizes is shown in Fig- 
ure ^. The PM-generated distributions agree at low flux 
values (high optical depths), but diverge for fluxes above 
0.8, with differences in the cumulative distributions as large 
as 0.05. As the resolution of the 9.6 Mpc box is in- 
creased to {Np,Ng) = (256,512), the amount of fine-scale 
structure continues to increase. The HPM-generated flux 
distributions differ markedly from the PM results, showing 
a smaller amount of fine-scale structure. The correspond- 
ing flux power spectra Pf (fc) are show in Figure ^ Again, 
as the grid resolution is increased in the 9.6 Mpc box, 
the amount of large-scale power continues to increase, while 
there is a decrease in the amount of power at shorter scales 
(fe > 0.04/iMpc"^ ). The 9.6 ft"^ Mpc HPM simulations 
show a deficit in power compared with the 9.6 Mpc PM 
simulations for k > 0.08/iMpc~^ , suggesting that the gas 
pressure forces have suppressed power at short scales. 
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Figure 2. Comparison of cumulative distributions of flux per 
pixel for the ACDM simulations for three different resolu- 
tions and three box sizes. The notation for each model is 
(Np , Ng ; boxsize; method), where the number of dark matter par- 
ticles is A'^, the number of grid zones is N^, and the boxsize is 
in units of Mpc. For the HPM simulations, the number of 
gas and dark matter particles are each A'^. The distributions di- 
verge at high flux values (small optical depths), suggesting small 
structures have not been resolved at the lower resolutions. 
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Figure 3. Comparison of flux Ppik) for the ACDM simulations 
for three grid resolutions and three box sizes. The notation for 
each model is as in Figure 2. A higher amount of power is pro- 
duced in the smallest box simulation compared with the others. 
The HPM model suppresses power at high k. 



Figure 4. Comparison of wavelet coefficient distributions for the 
ACDM simulations for 3 resolutions and 3 box sizes. 

3.2 Wavelet decomposition 

The cumulative distributions of wavelet coefficients are 
shown in Figure ^. The groups of curve from left to right 
correspond to the respective velocity scales 4 — 8kms^^ , 
8 - 16kms"^ , 16 - 32kms"^ , 32 - 64kms"^ , and 64 - 
128 km s"'^ . The small coefficients at the lowest scales show 
that no features are present on these scales. The physi- 
cally most interesting scales are represented by the 16 — 
32 km s~^ curves, corresponding to typical velocity widths 
of the absorption lines. The largest differences between the 
simulations are found on these scales. The convergence for 
the 9.6 Mpc box simulations is non-monotonic with 
and Ng, with the {Np,Ng) = (128,256) and (256,512) sim- 
ulations agreeing well, but not with the (256, 256) simula- 
tion. The latter however agrees closely with the (256, 512) 
12.48 Mpc simulation, suggesting that convergence has 
been reached with box size. The largest coefficients are found 
for the smallest box (6.24 fe"'^ Mpc ) simulation, indicating 
a greater amount of structure in the spectra on these ve- 
locity scales than found for the larger boxes. By contrast, 
the HPM simulation shows the least amount of structure on 
these scales. 



3.3 Line parameter decomposition 

A comparison of the cumulative distributions of A^hi and 
b determined by AutoVP is shown in Figure ^ for PM 
simulations for 3 different pairs of Np and Ng and 3 box 
sizes. As was found for the wavelet coefficient distributions, 
the behaviours for both the A^hi and b distributions are 
non-monotonic. The iNp,Ng) = (128,256) and (256,512) 
9.6h~^ Mpc simulations agree well in the b-distribution, 
but disagree with the {Np,Ng) — (256,256) simulations. It 
would appear that AutoVP is picking up unreal features 
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Figure 5. Comparison of cumulative distributions of absorption 
lino parameters for the ACDM simulations for 3 resolutions and 
3 box sizes. 



resulting from noise in the Np = 128 simulations that de- 
crease in going to Np = 256. Increasing the grid resolution, 
holding Np — 256 fixed, however, appears to recover real 
features that correspond statistically to the previously unre- 
solved fluctuations in the {Np,Ng) — (128, 256) simulations. 
Since the same trend is found for the wavelet coefficients, 
it seems unlikely that the non-monotonic behaviour is an 
artefact of the line finding and fitting procedures used by 
AutoVP. There are no obvious systematic trends with box 
size, the results otherwise appearing to have converged to 
better than 5%. 

Also shown are the resuhs for HPM for the 9.6 h'^ Mpc 
box and {Np,Ng) = (256,512). Both the column density 
distribution and the Doppler parameter distribution agree 
weU with the Np = Ng = 256 9.6 /i"^ Mpc PM simula- 
tion, demonstrating that incorporating the pseudo-pressure 
forces in the HPM simulation has very little effect on the 
resulting absorption line parameter distributions for this 
model. We note that the HPM simulation tends to pro- 
duce distributions shifted somewhat to higher column den- 
sity and Doppler parameter compared with the analogous 
{Np,Ng) = (256,512) PM simulation. 

A comparison of the results for HPM for different box 
sizes and particle/grid numbers is shown in Figure |^. The 
A^Hi distributions are well-converged. The fe-distributions 
would have appeared to converge, except the distributions 
shift toward larger 6-values in the highest resolution, small 
box simulation. Curiously, this is opposite the trends re- 
ported in the resolution studies of Theuns et al. and Bryan 
et al. , who found that the fe- values tended to decrease both 
with decreasing box size and increasing resolution, although 
at lower resolutions compared with our simulations. The 
trends found in the full hydrodynamics simulations are phys- 
ically sensible: a larger box adds longer wavelength power 
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Figure 6. Comparison of cumulative distributions of absorption 
line parameters for the HPM ACDM simulations for 2 grid reso- 
lutions and 2 box sizes, showing good convergence. 



to the simulations, which may increase the peculiar veloc- 
ities of the dark matter and the gas, while higher particle 
and grid numbers permit higher resolution of density fluc- 
tuations and shocks, which may result in a narrowing of 
absorption features. Possibly the opposing trend found here 
is a consequence of not allowing for shocks. For a typical 
sound speed of ~ 10 km s"'^ and peculiar velocities in the 
gas of ~ lOOkms"^ , shocks may be expected to play an 
important role in the velocity structure of the lines. 



4 COMPARISON WITH HYDRODYNAMICAL 
SIMULATIONS 

In this section we compare the results of the PM and 
HPM simulations against the Kronus simulation results of 
Machacek et al. (2000) (MOO) for a variety of cosmological 
models. The simulations were run for the same cosmological 
models as in MOO, matching the box sizes and the numbers 
of hydrodynamical cells and grid points for the gravitational 
force calculations. 



4.1 Flux distribution 

The distributions of flux for the 4 cosmological models 
considered are shown in Figure ^. Excellent agreement is 
found between HPM and the full hydrodynamical simula- 
tions (MOO), with the agreement least good for SCDM. The 
agreement with PM is not as good, but the cumulative dis- 
tributions still agree to within 10%. The agreement between 
both sets of A^p and Ng pairs for the HPM simulations shows 
that the results have converged at the relevant box sizes. 
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Figure 7. Comparison of cumulative distributions of flux per 
pixel for four cosmological models corresponding to the box sizes 
in MOO. The particle and grid numbers are shown as {Np,Ng). 

4.2 Wavelet decomposition 

The cumulative distributions for the wavelet coefficients are 
shown in Figure |^. Much better agreement is generally found 
between the PM simulations and the full hydrodynamical 
simulations than is found for HPM. The agreement is even 
poorer for the HPM simulations with {Np, Ng) = (256, 512) 
in all cases. The smaller values for the coefficients in the 
HPM spectrum indicates that the artificial hydrodynamics 
in HPM suppresses the full amount of velocity structure in 
the spectra compared with that found in the full hydrody- 
namics calculations. 



4.3 Line parameter decomposition 

The results of the Voigt line parameter analyses using Au- 
toVP are shown in Figures | to Both the PM and HPM 
results agree reasonably well with the results of MOO for 
the H I column density distributions. The best agreement 
for the Doppler parameter distributions is found for the PM 
simulations: the HPM simulations produce features that are 
too broad. While the PM and HPM simulations agree with 
the results of MOO to 5 — 10% on the TVhi distributions, the 
agreement between the Doppler parameter distributions is 
worse: ~ 10% for PM and ~ 20% for HPM. The disagree- 
ment between the Doppler parameter distributions tends to 
increase with the magnitude of the power spectra a,t z — 3. 



5 SUMMARY 

We have performed several simulations of the Lya forest 
using different background cosmological models, numerical 
codes and grid resolutions. Convergence to the distributions 



Figure 8. Comparison of wavelet coefficient distributions for PM, 
HPM, and the simulations of MOO (sohd lines). The PM results 
are shown for Np = Ng = 256 (dotted lines) and the HPM results 
for Np = 128, Ng = 256 (dashed lines). The groups of coefficient 
distributions are ordered as in Figure 4. 



of flux, wavelet coefficients, and absorption line parameters 
tends to be non-monotonic for the PM simulations as the 
number of particles and grid zones is increased, although 
agreements in the cumulative distributions to within 10% 
are achieved. 

A comparison with full hydrodynamical calculations 
shows that there is little advantage in HPM over PM, ex- 
cept possibly for the flux distributions. The HPM simula- 
tions generally result in broader absorption lines than PM. 
Both the PM and HPM simulations tend to yield broader 
lines than found in the full hydrodynamics calculations. This 
may be a consequence of the neglect of shocks in the PM and 
HPM simulations, which may sharpen the velocity structure 
of the absorption lines. 

The simulations show that the principal spectral prop- 
erties of the Lya forest may be qualitatively reproduced 
from the dark matter structures alone, to a remarkable accu- 
racy of 10% in the cumulative distributions of flux, wavelet 
coefficients, H i column density, and Doppler parameter. 
Tests of the models, however, require agreement at the 5% 
level per spectrum or better at the resolution of the Keck 
HIRES (Meiksin et al. 2000), so that ultimately full hydro- 
dynamical calculations will be required to confront models 
for the structure and physical properties of the Lyo? forest 
with direct observations. In order to match the observed 
spectra, it may ultimately be necessary to incorporate ra- 
diative transfer into the simulations both to allow for spatial 
fluctuations in the UV ionizing background and to reproduce 
the correct gas temperatures (and line widths), particularly 
if helium is completely photoionized at only moderate red- 
shifts. Several other physical effects may also be important, 
such as heating sources in addition to QSOs or galactic feed- 
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Figure 9. Comparison of line parameter distributions for PM, 
HPM, and the simulations of MOO for ACDM . 



Figure 11. Comparison of line parameter distributions for PM, 
HPM, and the simulations of MOO for SCDM. 
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Figure 10. Comparison of line parameter distributions for PM, 
HPM, and the simulations of MOO for OCDM. 



Figure 12. Comparison of line parameter distributions for PM, 
HPM, and the simulations of MOO for tCDM. 



back like winds. These effects may have sufficient impact on 
the statistics of the Lya forest that 10% accuracy is adequate 
to constrain them. Our results suggest that PM simulations 
alone may provide a reliable tool for investigating some of 
these issues without the considerable overhead of solving the 
full set of hydrodynamics equations. 
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Model 


Ho 






h 


n 


^8^-1 


ACDM 


0.4 


0.6 


0.0355 


0.65 


1 


1.0 


OCDM 


0.4 





0.0355 


0.65 


1 


1.0 


SCDM 


1 





0.06 


0.5 


1 


0.7 


tCDM 


1 





0.07 


0.6 


0.81 


0.5 



Table 1. Parameters for the cosmological models. Qq is 

the total mass density parameter, Qa the cosmological con- 
stant density parameter, Qf, the baryonic mass fraction, h = 
i^o/100 km s^-*^ Mpc^^ , where Ho is the Hubble constant at 
z = 0, n the slope of the primordial density perturbation power 
spectrum, and Cs^-i the fluctuation normalization in a sphere of 
radius 8h~^ Mpc. 

tional Center for Supercomputing Applications, University 
of Illinois, UrbanarChampaign. 



Model L(h-iMpc) A2(27r/L) Np Ng 



ACDM 


6.24 


0.64 


256 


256 


ACDM 


9.6 


0.45 


128 


256 


ACDM 


9.6 


0.45 


256 


256 


ACDM 


9.6 


0.45 


256 


512 


ACDM 


12.48 


0.35 


256 


512 


OCDM 


6.24 


0.98 


256 


256 


SCDM 


4.8 


0.50 


256 


256 


tCDM 


5.76 


0.19 


256 


256 



Table 2. Physical parameters for the PM simulations. L is the 

comoving box size in units of h^^ Mpc. Also shown is the dimen- 
sionlcss linear power spectrum A^(27r/L) at z = 3. The number 
of dark matter particles used is Np and the number of grid zones 
for the force calculation is Ng. 
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